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Abstract. We demonstrate the equivalence between the two popular models of nonlinear diffusion, the porous 
medium equation and the p-Laplacian equation. The equivalence is shown at the level of sclfsimilar solutions. 



1. Introduction 

The theory of nonhnear evolution equations of degenerate paraboUc type has been intensively studied in the last 
decades both for its importance in a number of applications in mechanics and heat propagation and also because it 
combines the theory of nonlinear evolution PDEs with geometry in the form of free boundaries. Two of the most 
studied models have been the porous medium equation, shortly PME, and the p-Laplacian equation, shortly PLE. 
We write the PME as 

(1) ut = A(u™/m), m > 1. 

It generalizes the heat equation, which is the linear case m — I, and we also include the range of parameters m < 1 
where it is known as the fast diffusion equation, [TU] (m < is also allowed in the general theory). On the other 
hand, the standard PLE is 

(2) Ut = ApM := V • (I V^r^^ V"), 

where the parameter ranges in the interval 1 < p < oo. The heat equation is the case p ~ 2. There has been strong 
interest in recent times in the case p = 1 as a geometrical flow, [1] . 

It was soon remarked that the theory of both equations offers striking parallels. To mention just one, the PME 
has the property of finite speed of propagation for m > 1, while the PLE has the same property for p > 2. The 
consequence in both cases is that solutions with compactly supported initial data u(x, 0) > preserve the property 
of compact support w.r.t the space variable for all t > 0, and this implies the existence of a clear-cut interface or 
free boundary separating the non-empty regions where m > and where u = 0. Many other properties are known 
that reinforce that parallelism, cf. [5] , [5] . 

In this note we show that the very strong connection between the two equations can be improved into complete 
equivalence when we consider special classes of solutions. Such solutions are the self-similar solutions, i.e., solutions 
of one of the two forms, 

(3) uix, t) = t-^f{x r^) = t-'^fir,), u{x, t) = {T- tr,f{x [T - tf) = (T - 

We will call the solutions of the form Q-left sclfsimilar solutions of Type I, those of the form Q-right are called 
Type II. The practical importance of such solutions lies in the fact that they usually describe the asymptotic 
behavior of solutions of Cauchy problem associated to parabolic equations, as the extensive literature shows, cf. 
[3], [TU] and their references. In view of the definite progress obtained recently in the study of the PME, as shown 
in the latter reference, and the key role actually played by suitable sclfsimilar solutions in the classification of the 
basic qualitative and asymptotic properties, the equivalence is of clear significance to establish a similar theory for 
the PLE. We remark that for our equations there is still a possible selfsimilarity of Type III, with solutions of the 
form u{x,t) — e"*/(xe'^*), but it plays a minor role in the theory. Since the details are almost the same we will 
make only side mention to it. 

A popular technique used for the description of these solutions is phase-plane analysis, cf. [3l O [1^. Our 
main contribution in this Note consists in introducing a new set of variables in the phase-plane description of the 
self-similar solutions of the p-Laplacian equation. They are not straightforward at all, but they have the property 
that the resulting phase plane can be exactly identified with a known version of the phase-plane for the self-similar 
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solutions of the PME, the only apparent difference amounting to different formulas for the values of the parameters. 
Here is our main result. 

Theorem 1.1. The analysis of radial selfsimilar solutions for both the PME (when m ^ 1) and the PLE (when 
p ^ 2) can be reduced to a particular case of the autonomous ODE system 



(4) 



$ ci$^ - C2*<I> - C3$ ± * + sgn(b) 



where the parameters ci, C2, C3 and b are explicit functions of n, m and (3 in the PME case; of n, p and (3 in the 
PLE case. The variables $ and ^' are given by algebraic expressions in terms ofri,f,f', which are different for 
both equations. It follows that a translation rule can be set up from the PME into the PLE and vice versa, which 
involves changing the space dimension in the process. Besides, there are in general two options for the equivalence 
map in either direction. 

Obviously, in the exceptional case m ~ 1, p ~ 2 the two equations coincide and the identity transformation 
solves the equivalence problem. We recall that there exist a large number of selfsimilar analyses of the PME and 
the PLE but the equivalence has not been remarked. 

In the sequel we give whole details of these assertions: a careful selection of the correspondence between the two 
sets of parameters allows for a complete identification of the phase planes, which implies a transformation rule for 
the set of selfsimilar solutions from one equation to the other. We then derive some interesting applications. We 
will also provide self-maps of the solutions of the PME (resp. the PLE) into themselves based on an repeated use 
of the transformations, using a different option at each turn. We point out that when performing the selfsimilar 
analysis with radial profiles by means of ordinary differential equations we may assume that the space dimension 
is any real positive number. 

2. Phase plane analysis for the PME 

We want to find solutions in the self-similar forms ^ for equation ([T]) posed for x G with m 7^ 1. This is 
well-known, and is described in detail in [9] for instance. First, the relation between the similarity exponents a 
and P reads (m — l)a + 2/3=1 for Type I selfsimilarity, and (m — l)a + 2(3 = —1 if it belongs to Type IL Then, 
under the usual assumption of radial symmetry the profile / must satisfy the ODE 

(5) ryi-"(r;"- V"- V')' + «/ + f^vf = 0, 

where rj > 0. From now on we assume that n is any positive real number. We first introduce the variables: 

(6) X = r^r/f, Y = r^^f'-^. 

We also replace the 77 variable by r ~ \ogi]. The functions X{r) and Y{r) satisfy the autonomous ODE system: 

(7) X = (2 - n)X - mX^ - (a + f3X)Y, F = (2 + (1 - m)X)Y. 

where the over-dot indicates differentiation with respect to r. Now, we introduce a new pair of variables 

$ = (2 + (1 - m)X)/v4^, ^ = Y/\b\, 

where b ~ 2n{m — m^jijfi — 1) and 777c = ~ 2) jn. This assumes that 777 7^ 777c so that & 7^ 0. Replacing then the 
r variable by r\ ~ \/[6[r, so that over-dot indicates differentiation with respect to ri, the system takes the desired 
quadratic form with precise values for the constants given by 

m a /Tu {n + 2){m,- ms) n-2 

m-1 (^rn -l)y^\ n + 2 

and the ± of the equation is -I- for the Type I and it is — for the Type II. With these values System (jlj has free 
parameters m, n and (3, since a can be calculated from them. The case ?77 = ?77c will be discussed below. 
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3. Phase plane analysis for the PLE 

In this case, the relation between the similarity exponents a and (3 is {p — 2)a+pf3 = 1 for Type I selfsimilarity, 
and {p — 2)a + pf3 = —1 if it is of Type II. The radially symmetric profile / satisfies the ODE 

(8) v'-'^W'-'lfT'iy + af + IM' = 0, 

where prime denotes differentiation with respect to > 0. From now on we assume that n is any positive real 
number. In a similar way as before, for p ^2 we introduce phase plane variables, a bit different from the ones in 
the PME case: 

(9) X^-r/^l/'ji-^'/', Z = i-P!, where 7=^7^ 

2 — p 

We thus get the autonomous ODE system 

(10) ^^-^X^-(n--i)X + aZ\X\-^ -^\X\X, Z = jZ -\X\^X, 
2-p 

where we have replaced the rj variable by r = logr/, so that over-dot indicates differentiation with respect to 
r. This system is not quadratic so that we perform further change with this objective in mind. We introduce 
Y = \X\p^ XZ = —i]\f'\^^f'f ^-nd the flow equations become 

(11) X ^'^—^X h-n + aY ~ PIXI) , Y = -aY^ + nY + I3Y\X\ - \X\, 

p-1 

which is a quadratic system if X has a sign. For the next step we assume that X > and we set 

^ = aX, $ = ^ '^^ (j-n + aY - BX) , 

ip-l)^m 

where 

1 p{n+l){p-pc) 2n 
"= ~r :^T7 7^ Pc = 



|6|(p-l)' {p~2){p-l) n+1- 

In order to proceed further we assume that p ^ pc- After all these transformations, the flow equations become 
exactly the desired The ± of the equation is + for the Type I and it is — for the Type 11, and the constants 
are now given by 

Cl = C2=P\/\b\, C3 = - - 



Ps 



p-2 ip-2)^\ ' n + 2 

This system has free parameters p, n and /3, and p = 2 is excluded. We have replaced the r variable by ri = A/j^fr, 
so that over-dot indicates differentiation with respect to ri . 

The critical cases, tUc and pc- Some changes have to be made in the special cases m = rric for the PME or 
p = Pc for the PLE, since our definitions imply that 6 = 0. we change it into y/b = (n — 2) for PME and \/b = n for 
PLE so that C3 = — 1 and the independent term sgn{h) disappears from the second equation of system ([4]), which 
becomes 

(12) ^' = ^f^, $ = Ci$^-C2*$ + $±^'. 

with Cl and C2 as before. 

Tpe III selfsimilarity. It offers few novelties. The exponent relationships are a{l — m) ~ 2(3 and a{2 — p) ^ pf3 
respectively. We obtain a similar system, except for the fact that the term ±5* disappears. Summing up, the 
coefficient of is + for Type I, — for Type II, and for Type III. 

4. Equivalence 

Identifying the two systems implies taking parameters (m, n, 0) in the first system and (p, n' , /?') in the second, 
so that the expressions for the free constants of System (jl]) obtain the same values, (i) Identification of ci gives 
the relation p = to + 1, which is well-known from different dimensional considerations, (ii) Identifying ca gives 

{n + 2){m-m,) {n' + 2){p-p,) 



(to-I)VH (P-2)V|6'| 
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Note that b' = p{n' + l){p- Pc)/{p - 2){p - 1) with pc = 2n'/{n' + 1). This leads to 

2m{n{m - 1) + 2){n{m - 1) + 2{m + l)f = (m + l)(n'(m - 1) + to + l)(n(m - 1) + 2{m + 1))^, 
which can be written as a quadratic condition: ^(tt,')^ + Bn' + C = 0, with 

A = 2m(m - \ f{n{m - 1) + 2], B = -(to + 1)(to - \f{n - 2f, C = -(m - lf{m + lf{n - 2f . 

This gives the values for n' as a double function of n, representing the change of dimension from PME to PLE 

, {n-2){m + l) , {n-2){m+l) TOc(to+1) 

(13) = 7, . "2 = 7, = • 

Zm n — 2 — nm rUc — m 

The double solution exists for m ^ 0,mc or 1. The two values of n' coincide for to = to^, and then n' = n and 
p — Ps- Clearly, we have 

1 I _ B_l-m_2-p 
n[ n'2 C 771+1 p 
which expresses a possible self-map of the PLE (for fixed p) by change of dimension. The situation is similar in 
the other sense: if p 7^ Pc{n') = 2n'/(n' + 1), there are two possible values of n for every given n' and we have the 
relationship 

1 1 1 - TO 



rti — 2 712 — 2 2to 

which expresses a possible self-map of the PME. A precedent of self-relations can be found in King's [5]. 
iii) Identifying C2 implies that /3-\/[6[ = 13' ^J\h'\, so that 

^ ^ {m-l)p{n' + l){p-p,) 
(3"^ 2r7,(TO — TOc)(p — 2)(p — 1) 
which gives the corresponding two values for /3' in terms of (3 

l3[ 2m n-2 n(m - 1) + 2 2-77 



/3 TO + 1 n'l (3 TO + 1 772 

Note that in both cases I3'^{n - 2)^ = {[3'n'Y. 

(iv) We still have to check that sgn{b) ~ sgn{b'). We have 

p{n' + l){p- pc){m- 1) ((to - 1)77' TO+ 1)(to + 1) (77')^ 



b'/b 



2n{p - 2){p - l){m - rric) 2m{n{m - 1) + 2) (77-2)2' 

We have used that p = to +1, pc ~ 2n' / (77' -1- 1), and (to — 1)77' -f 777 + 1 = 2m{n{m — 1) -f 2){n'Y / (to +\){n — 2)^ 
given by the quadratic condition for 77'. 

Some general conclusions. (1) The transformation is not applicable for 77 = 2 since it implies 77' = for both 
branches. Besides, the study is different for the cases 77 < 2 and n > 2. 

(2) In the case 77 = 1 we have same dimension 712 = 1 for all to and then (3' = j3. This is well known since formally 
differentiating the one-dimensional PLE gives the PME. For the other branch corresponding to n = 1 we have 

TO + 1 (3 



77i 



2^ r' / 7 

TO 77]^ 

so that 77^ > if — 1 < TO < 0, a case of very fast diffusion studied in [4], see also [9]. Also, 77'^ > 1 if —1/3 < to < 0. 

(3) Cases 77 > 3. (i) We have n' = 77 if and only if to = nis in both branches. This is the Yamabe case and 
it allows us to get an explicit equivalent of the Loewner-Nirenberg solution of the fast Diffusion Equation for the 
p-Laplacian equation. 

(ii) In the case of the first branch, 77' > if to > or to < —1 (a case we may disregard since then p < 0). In the 
first branch 77' decreases with m from 00 as to — > to (?7 — 2)/2 as m — > 00. 

(iii) In the second branch n' is positive for — 1 < to < TOc = (77 — 2)/77 and increases from to 00 in that interval; 
more precisely, 77' = 1 for to 

(4) Identification of the critical cases 777 = 777c in the PME and p = pc in the PLE takes place for 77' = 77 — 1 and 
/3' = p{n - 2)/n' = [3{n - 2)/(77 - 1). 

(5) In the limit 777 ^ 1, p ^ 2, both equations reduce to the linear heat equation (there is also the possibility 
of converging to the eikonal equation ut ~ jViip with a different limit process, cf. [5]). Now, extrapolating our 
transformation to the linear case, to = 1, p = 2, we have n'l = n — 2 with (3' = [3, and 773 = 2 — 77 with (3' = (3. 
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5. Application. Some explicit solutions 

Let us now analyze some self-similar solutions given by the plane of phases for the equation of porous medium, 
and compare these solutions with the corresponding solution of the p-Laplacian equation. We will be mainly 
interested in dimensions n > 3 but lower dimensions appear in the final results. 

Barenblatt solution and straight lines. Consider the equation of the trajectories with sgn{h) — 1 and 
selfsimilarity of Type I: 

d$ _ Ci$^ - 02*$ - C3$ + + 1 

^ ' rf* ~ ^$ ■ 

We want to find the linear solutions to (I14p . In order to do that we assume that $ = ai^" + 02 and insert this form 
into (jl4p . obtaining the following conditions: ai = 02 = C2/(ci — 1), and C2 must be a solution of the equation 
C1C2 — (ci — l)c3C2 + (ci — 1)^ = 0. Since C2 = P\/W\^ two values for (3 and (3' which are 

(15) (3i = — ^ or /32 = 77- (porous medium) 

n{m — Ij + 2 Am 

(16) P'l = —r, \n , or /32 = - (p-Laplacian) . 

The profiles corresponding to the trajectories for /?i and (i'^ produce the so-called Barenblatt solution, both for the 
PME and the PLE. The other case is as follows: (i^ corresponds to the dipole solution of the porous medium, 

1 

/ 1 \ 

f = ■q^~^ Ik ^(mn-n+2)/va 

ci. [To], Section 4.6, and /Jj corresponds a'2 = 0, and a profile whose derivative is 

l/(p-2) 



+ 



Yamabe case. It corresponds to ps = rus -I- 1 = 2n/{n + 2) (cf. [9], Section 7.5 and Appendix AIII). Choosing 
Type II and putting /3 = 0, so that the self-similar solution is a separate variables solution, we have ci = — (n — 2)/4, 
C2 = C3 = and sgn{h) = 1. The equation of the trajectories turns out to be 

d$ _ ci*^ - * + 1 

Solving the above Bernoulli's equation, we obtain that the curve "if = n(l/(n — 2) — 1/4$^) is one of its solutions. 
The explicit profiles that we get after undoing all the changes of variables are: 

(17) f(ri) = iki -f - — — ) 2 (Loewner-Nirenberg profile for the PME with m = m^) 

^ _„ 4n /4n3fc2 A'"'^^^"* 

(18) fM ~ C(l + k2'r]"-^ ) 2 ^ C = ^5 I (new p-Laplacian profile for p — ps) 

n + 2 — A J 

where ki , fc2 > are arbitrary. 

Note. A detailed account of the results and some applications will appear as a separate publication [7]. 

References 

[1] A. Andreu-Vaillo, V. Caselles, J. M. Mazon. Parabolic quasilinear equations minimizing linear growth junctionals, Birkhauser, 
Basel, 2004. 

[2] D. G. Aronson, J. L. Vazquez. The porous medium equation as a finite-speed approximation to a Hamilton-Jacobi equation, Ann. 

Inst. H. Poincare Anal. Non Lineaire 4 (1987), no. 3, 203-230. 
[3] G. L Barenblatt. Scaling, Self-Similarity, and Intermediate Asymptotics, Cambridge Univ. Press, Cambridge, 1996. Updated 

version of Similarity, Self-Similarity, and Intermediate Asymptotics, Consultants Bureau, New York, 1979. 
[4] J. R. Esteban, A. Rodriguez, J. L. Vazquez. A nonlinear heat equation with singular diffusivity. Comm. Partial Diff. Eqs. 13 

(1988), 985-1039. 

[5] M. F. Bidaut-Vcron. Self-similar solutions of the p-Laplace heat equation: the fast diffusion case. Pacific Journal of Mathematics, 
227 (2006), 201-270. 

[6] E. DiBenedetto. Degenerate Parabolic Equations, Series Universitext, Springer- Verlag, New York, 1993. 
[7] R. lagar, A. Sanchez, J. L. Vazquez, paper in preparation. 



6 



ARIEL sAnCHEZ AND JUAN-LUIS VAZQUEZ 



[8] J. R. King. Self-similar behaviour for the equation of fast nonlinear diffusion, Phil. Trans. Roy. Soc. London A 343 (1993), 
337-375. 

[9] J. L. Vazquez. Smoothing and Decay Estimates for Nonlinear Diffusion Equations. Equations of Porous Medium Type. Oxford 
Lecture Series in Mathematics and Its AppUcations 33, Oxford University Press, 2006. 
[10] J. L. Vazquez. The Porous Medium Equation. Mathematical Theory, Oxford Mathematical Monographs, Oxford Univ. Press, 
2007. 

E-mail address: ariel . sanchezOur jc. es 
E-mail address: juanluis . vazquezauam. es 

A. S.: Departamento de Matematica Aplicada, Universidad Rey Juan Carlos, Mostoles. 28933 Madrid, Spain 
J. L. v.: Departamento de Matematicas, Universidad Autonoma de Madrid, Cantoblanco. 28049 Madrid, Spain. 



